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Entanglement is the resource to overcome the natural limitations of spatially separated parties 
restricted to Local Operations assisted by Classical Communications (LOCC). Recently two new 
classes of operational entanglement measures, the source and the accessible entanglement, for ar¬ 
bitrary multipartite states have been introduced. Whereas the source entanglement measures from 
how many states the state of interest can be obtained via LOCC, the accessible entanglement 
measures how many states can be reached via LOCC from the state at hand. We consider here 
pure bipartite as well as multipartite states and derive explicit formulae for the source entangle¬ 
ment. Moreover, we obtain explicit formulae for a whole class of source entanglement measures that 
characterize the simplicity of generating a given bipartite pure state via LOCC. Furthermore, we 
show how the accessible entanglement can be computed numerically. For generic four-qubit states 
we first derive the necessary and sufficient conditions for the existence of LOCC transformations 
among these states and then derive explicit formulae for their accessible and source entanglement. 


I. INTRODUCTION 

Entanglement is considered to be one of the charac¬ 
teristic traits of quantum mechanics. Besides its inter¬ 
est from the foundational point of view, it plays a key 
role in quantum information science, being a resource 
for most of its applications such as quantum communica¬ 
tion and quantum computation [1]. This has led to the 
development of entanglement theory [2], which has also 
recently brought new tools for the understanding of quan¬ 
tum many-body systems in condensed matter physics [3] . 
The characterization (of different forms) of entanglement 
and its quantification play a central role in this theory. 
Ideally, entanglement measures should allow to opera¬ 
tionally quantify how useful a state is for quantum infor¬ 
mation applications, how efficient protocols can be and 
provide quantitative means to link entanglement to other 
physical phenomena such as phase transitions. Despite 
remarkable advances, there are still many issues that re¬ 
quire a better understanding, particularly in the case of 
multipartite states and mixed states in general. There 
is a clear need for entanglement measures which are on 
the one hand operationally meaningful and on the other 
hand easy to compute. 

The paradigm of local operations assisted by classical 
communication (LOCC) is of paramount importance in 
the theory of entanglement. LOCC constitutes the most 
general class of transformations allowed by the rules of 
quantum mechanics to distant parties: each subsystem 
can undergo any form of local quantum dynamics (i.e. a 
completely positive map) at will and the parties’ actions 
can be correlated through the use of classical communi¬ 
cation. Thus, LOCC maps constitute all possible proto¬ 
cols for the manipulation of entangled states. Moreover, 
LOCC allows to formulate entanglement theory as a re¬ 
source theory. This is because entanglement cannot be 
created by LOCC alone and it is, hence, a resource to 


overcome the limitations of parties restricted to this class 
of operations. The basic law of entanglement is that it 
cannot increase by LOCC and, therefore, an entangle¬ 
ment measure must be a function that is monotonic with 
respect to the ordering induced by this class of transfor¬ 
mations. Unfortunately, this ordering is in general just 
partial (see e.g. [4]) and many different functions that 
qualify as entanglement measures exist. For this reason, 
it is desirable to seek for some further operational mean¬ 
ing behind these functions in the context of LOCC in 
order to single out the most relevant entanglement mea¬ 
sures. 

The identification of the operational meaning of entan¬ 
glement measures has been particularly successful in the 
case of pure bipartite entangled states. It has been shown 
that the entropy of entanglement represents the rate at 
which any such state can be reversibly transformed by 
LOCC in the asymptotic limit of infinitely many copies 
into the maximally entangled state [5]. Here, the notion 
of maximal entanglement stems from LOCC transforma¬ 
tions on single copies: the maximally entangled bipar¬ 
tite state is the unique state that cannot be obtained 
from any other Local Unitary (LU)-inequivalent bipar¬ 
tite state but allows to be transformed into any other 
state [4]. Regrettably, the idea of reversible entangle¬ 
ment distillation cannot be extended to bipartite mixed 
states. The entanglement cost (i.e. the rate at which the 
state is obtained from the maximally entangled state) 
is in general different from the distillable entanglement 
(i.e. the rate at which the state is transformed to the 
maximally entangled state) [6]. Furthermore, the com¬ 
putation of these measures is formidably difficult. There 
are other measures that boil down to the entanglement 
entropy for pure states but which are different for mixed 
states such as the relative entropy of entanglement [7] or 
the squashed entanglement [8]. However, despite inter¬ 
esting properties, they lack the interpretation in terms 
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of LOCC-conversion rates and their exact computation 
is still very hard. 

The situation is even worse in the multipartite case. 
There exists no unique maximally entangled state [9] and 
it is not clear to which set of states LOCC asymptotic 
rates have to be defined (the so-called MREGS (minimal 
reversible entanglement generating set) problem [10]). 
Nevertheless, several entanglement measures have been 
proposed from the purely mathematical perspective of 
invariants such as the tangle [11] to the ability to create 
bipartite entanglement among given cuts such as the lo- 
calizable entanglement [12]. Besides this, many bipartite 
measures admit a generalization to multipartite states 
like distance-based measures such as the geometric mea¬ 
sure of entanglement [13] to cite an example. However, 
all these measures are hard to compute and, more im¬ 
portantly, many of them do not have an interpretation 
in the context of the LOCC paradigm. 

In order to try to close this gap, we have recently in¬ 
troduced in [14] operational entanglement measures for 
general (pure or mixed) multipartite states which are 
LOCC-meaningful. The basic idea is to focus on the ca¬ 
pacity of a state for single-copy LOCC transformations 
instead of transformation rates in the asymptotic regime. 
In more detail, we have introduced two classes of entan¬ 
glement measures based on the accessible volume and the 
source volume of a state. The first one quantifies the rel¬ 
ative volume of inequivalent states that can be accessed 
from our state by LOCC while the second one quantifies 
the relative volume of inequivalent states from which our 
state can be obtained by LOCC. Hence, these measures 
have a clear operational meaning in the context of LOCC 
transformations. The larger the accessible volume is the 
more entangled (i.e. useful) a state should be: the state 
is at least as powerful in applications as any state in the 
accessible set as any protocol achievable with the latter 
can also be achieved with the former by converting it by 
LOCC to this accessible state first. On the other hand, 
the larger the source volume the less entangled a state 
should be given that many inequivalent states can reach 
this state by LOCC. Notice then that these measures are 
completely general, valid for arbitrary states of any di¬ 
mension. Moreover, they can be generalized to classes of 
entanglement measures by considering different Hilbert 
spaces for the initial and the final state (see [14] and Sec. 
H). Furthermore, whenever the possible LOCC transfor¬ 
mations are characterized, these measures can be com¬ 
puted. Although characterizing LOCC transformations 
is notoriously difficult [15], this problem is experiencing 
some progress lately [4, 9, 16]. Using these insights, in 
[14] we have analyzed the 3-qubit pure-state case provid¬ 
ing explicit formulae for these measures. Here, we fur¬ 
ther show the versatility of this approach by computing 
these measures for the bipartite pure-state case and the 
generic four-qubit pure-state case. In order to do so, we 
complete the analysis of [9] and characterize all possible 


LOCC transformations among the latter class of states. 

The outline of the remainder of the paper is the fol¬ 
lowing. In Sec. H we recall the definition and properties 
of the two classes of entanglement measures we intro¬ 
duced in [14]. In Sec HI we consider general bipartite 
pure states and derive explicit formulae for the source 
entanglement and its generalizations. That is, we obtain 
explicit formulae for a whole class of operational entan¬ 
glement measures that quantify how easy it is to gener¬ 
ate a bipartite pure state from other bipartite quantum 
states via LOCC. These measures can be used to charac¬ 
terize, e.g., the entanglement contained in pure quantum 
states of two qubits or two qutrits. Moreover, we demon¬ 
strate how the accessible volume can be computed and 
illustrate the results by considering up to two four-level 
systems. Sec. IV deals with the multipartite case, where 
we first derive the necessary and sufficient conditions for 
LOCC convertibility of generic four-qubit states and then 
compute the new entanglement measures. 


II. ACCESSIBLE AND SOURCE 
ENTANGLEMENT 

In this section we review the definition of the entan¬ 
glement measures we have introduced in [14]. Since lo¬ 
cal unitary (LU) transformations are reversible LOCC 
transformations, the entanglement of two states is equiv¬ 
alent if they are related by such a transformation. Hence, 
we consider the possible LOCC transformations among 
the LU-equivalence classes of states rather than states in 
general. That is, we always pick a unique representative 
state from each class. We say that a state p can reach 
a state a and that a is accessible from p if there exists 
a deterministic LOCC protocol which transforms p into 
cr. For a given state, p, we denote by Ma{p) the set of 
states which can be accessed from p and by Ms{p) the 
set of states which can reach p (see Fig. 1). Let p, de¬ 
note an arbitrary measure in the set of LU equivalence 
classes of states. Then, the source volume of p is de¬ 
fined by Vs{p) = p[Mg{p)\ and the accessible volume by 
Va{p) = p[Ma{p)]. As mentioned in the introduction, 
these quantities measure respectively the amount of in¬ 
equivalent states that are not less useful and not more 
useful than the state at hand. Thus, the accessible en¬ 
tanglement and the source entanglement are defined by 

J 71 / _ ^(P) jp / 

^a\P) ^sup ; ^s\P) 1 

where denote the supremum of the accessible 

(source) volume according to the measure p. 
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FIG. 1: (color online). In this schematic figure the source 
set, Ma{p), and the accessible set, Ma{p), of the state p are 
depicted. Any state in Ms{p) can be transformed to p via 
LOGO and p can be transformed into any state in Ma{p) via 
LOCG. 

A few remarks are in order. First, due to their op¬ 
erational meaning, it is easy to see that Ea and Eg are 
valid entanglement measures (see [14]). That is, for any 
deterministic LOCC protocol A it holds that 

E,{A{p)) < Ea{p), Eg{A{p)) < Es{p) (2) 
for every state p. 

Notice that the so-called entanglement monotones have 
been sometimes considered as entanglement quantifiers. 
These quantities must, in contrast to entanglement mea¬ 
sures, fulfill that if different outcomes of the LOCC map 
A acting on a state, Afe(p), can be postselected, each of 
them occurring with probability pfc, then entanglement 
cannot increase on average, i.e. J2kPkEiAk{p)) < E{p). 
However, we follow here the definition of entanglement 
measures, for which E{A{p)) < E{p) for any LOCC map 
A must hold. For a discussion justifying the suitability 
of this requirement over the averaged one, we refer the 
reader to Sec. XV.B.l of [2]. 

Second, the choice of the mathematical measure p to 
compute the volumes of the source and accessible sets is 
in principle arbitrary. That is, any valid measure leads 
to an entanglement measure. Depending on the measure, 
one might obtain a different ordering of the states, as is 
common in case of larger systems than those composed of 
two qubits [17]. The choice of the measure can be fixed 
by the physical constraints of the problem at hand or 
even out of mathematical convenience. Obviously how¬ 
ever, one should always consider a fixed dimension of the 
states in the source set. This is because, by imposing no 
constraint on the dimension, the size of this set would 
explode by considering source states of arbitrarily large 
dimension. A similar thing happens for the accessible 
set as it would shrink to measure zero relative to states 
of sufficiently large dimension. Thus, one should always 
restrict to a choice of measure which is supported on 
states of a given dimensionality. Note that, by changing 
this choice of dimension (and/or the number of consid¬ 
ered subsystems), one has the freedom to obtain different 
families of measures {Eg}k>d and {E^}k<d, where d is 
the effective dimension of the subsystems of the state at 
hand [14]. 


Another issue to take into account is that it can be the 
case that one is interested in the relative volumes of states 
under LOCC transformations among particular classes of 
states. In this case one should choose measures which are 
only supported on these classes. For instance, and this 
will always be the case in this article, one can consider 
only transformations among pure states. Hence, the mea¬ 
sure p will only be supported on the set of LU-equivalence 
classes of pure states. This seems a reasonable choice as 
the study of pure entanglement inside the full (mixed) 
state space might be too coarse grained given that pure 
states are of measure zero (e.g. Mg would have measure 
zero for all pure two-qubit states on the full state space 
of two-qubit states given that it is not possible to convert 
via LOCC a non-pure two-qubit state into a pure entan¬ 
gled two-qubit state [18]). Another instance is the study 
of multipartite states where we will consider measures 
that are not supported on biseparable states. Finally, it 
might be even desirable to choose different mathematical 
measures for different subclasses. This might be the case 
when these subclasses are not LOCC related (thus, this 
does not compromise the validity of Eq. (2)) or when the 
different subclasses have source or accessible sets living 
in manifolds of different dimensionality (see [14]). 

III. BIPARTITE SYSTEMS 
A. Preliminaries 

1. LOCC transformations for bipartite systems 

Every pure state of a bipartite quantum system with 
Hilbert space % = ® can be (up to local uni- 

taries) written as ]■)/) ^lu where d = 

min{di,(i 2 } and > 0 denote the Schmidt coefficients 
with = 1- We denote by = (Ai,..., Xd) € R‘^ 

the Schmidt vector of {ip) and will consider in the fol¬ 
lowing w.l.o.g. two d-level systems. As can be easily 
seen from the Schmidt decomposition given above, two 
d-level states, J^), \cj)) , are LU equivalent if and only if 
(iff) X^{ip) = A4((/)), where here and in the following 
X^{'4)) € R‘^, with X^{'ip)i > A-*'(^/)i+i > 0 denotes the 
sorted Schmidt vector of j^). In the context of LOCC 
transformations of pure bipartite states, the following 
functions of a; = (xi,..., Xd) G 

k 

Ek{x) :='^Xi, k G {l,...,d}, (3) 

play an important role. It was shown in [4] that a state 
\tp) G fi can be transformed into \4>) G 'H determinis¬ 
tically via LOCC iff A('0) is majorized by A((/), written 
A(V') ^ X{(j)), i.e. 

EkiX^Pj)) < EkiX^cf)) Vfc e {1,..., 4, (4) 
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with equality for k = d. 

A direct consequence of this criterion is that the source 
and accessible set of |'0) GH are given by 

MsW = € n s.t. X{(j)) < (5) 

Maitp) = {|0) e y. s.t. A(V’) ^ X{4))}. (6) 

As mentioned above, the volumes of these sets measure 
how easy it is to generate a given state and how many 
states a given state can access respectively. As explained 
before, we quantify the set of LU-equivalence classes as 
the amount of LU-equivalent states which can be used to 
reach a given state is not of relevance here. Due to the 
one-to-one correspondence between the LU-equivalence 
classes of bipartite pure states and the sorted Schmidt 
vectors, we can associate to the sets given in Eqs. (5) 
and (6) the following sets of sorted Schmidt vectors in 
R^: 

MsW = {A^ G R‘^ s.t. A^ ^ A(V’)}, (7) 

MaW = {A^ G R‘^ s.t. A(V') ^ A^}. (8) 

Unless otherwise stated, d denotes the Schmidt rank of 
IV') and these sets are hence supported on states of the 
same dimensions as \tjj). Notice, however, that, as ex¬ 
plained in the previous section, it is possible to consider 
source sets supported on states of larger dimension and 
accessible sets supported on states of smaller dimension. 
This gives the possibility to obtain different families of 
measures (see [14] and Sec. IIIB below). 

The sets given in Eqs. (7)-(8) are convex polytopes. 
In order to compute their volumes in the subsequent sec¬ 
tions, we first review some of their properties. 

2. Convex polytopes 

We briefly recall some definitions and results concern¬ 
ing convex polytopes. The reader is referred to [19] for 
a comprehensive introduction into the study of convex 
poly topes. A closed halfspace in is a set of the form 
H = {x € R^ s.t. cx + Co > 0}, with c G 1 ?^, cq G R. We 
denote by h{H) the hyperplane that fulfills the inequality 
with equality, i.e. h{H) := {x G R^ s.t. ca; -|- cq = 0} . 
A well-known and important result for subsets P C R!^ 
is the following equivalence. A subset P C R^ is the 
bounded intersection of a finite set of closed halfspaces, 
i.e. 

P = {x a R^ s.t. Ax -I- 6 > 0} for A G R^^^, b G i?™, 

( 9 ) 

iff it is the convex hull of a finite point set, i.e. 

P = conv(U), where V = {vi G R’^}. (10) 


The object P is called a convex polytope and a repre¬ 
sentation of the form in Eq. (9) is an i7-representation 
of P. The latter way to define the polytope is called V- 
representation. It is unique if the set {ui} is minimal. 
This minimal set is called vertex set, vert(P), of P. The 
dimension of P, dim(P), is the dimension of its affine 
hull, i.e. of the smallest affine subspace of R^ that con¬ 
tains P. A face P of P is a set of the form P = Pr\h{H) 
for some hyperplane h{H). Here, H D P is a closed 
halfspace that contains P. The dimension of F is again 
the dimension of its affine hull. Vertices are faces of di¬ 
mension 0, while faces of dimension 1 and dim(P) — 1 
are called edges and facets of P C respectively. Fur¬ 
thermore, a polytope with dim(P) = A: is called simple 
if every vertex is contained in only k facets. Note that 
this is the minimal number of facets a vector must be 
an element of in order to be a vertex. That is, every 
vertex of a simple polytope is an element of only k hy¬ 
perplanes that are associated to closed halfspaces in a 
minimal P-representation, i.e. an P-representation of 
the polytope that consists of the minimal number of in¬ 
tersecting halfspaces. Two vertices are called neighbors 
if they are connected by an edge, i.e. if fc — 1 of the above 
mentioned hyperplanes coincide. It is easy to see that a 
/c-dimensional polytope is simple iff every vertex has k 
neighbors [29]. 

Whereas sets such as Ms{tp) have been studied in de¬ 
tail in the literature (see [19] and references therein), and 
its vertices have been computed [20] (see Sec. IIIB), 
much less is known about the accessible set. Note that 
not even the identification of the vertices of the accessi¬ 
ble volume is straightforward. In fact, even the number 
of vertices depends strongly on the state of interest, as 
we will explain in Sec. IIIC. However, there exist algo¬ 
rithms that can compute the vertices, i.e. the minimal 
U-representation, of a convex polytope from a given P- 
representation. They can also be used to compute the 
accessible volume, which we will use in Sec. HI C to de¬ 
termine the accessible entanglement. 

B. Source volume 

We determine here the volume of the source set given 
in Eq. (7). In order to do so, we give a description of 
AisO^) in terms of convex geometry. Let us first note 
that the conditions in Eq. (4) on the sorted Schmidt 
vectors are equivalent to the following set of inequalities 
on the unsorted Schmidt vectors 

n = {E^ixm < EkiXH'^)) s.t. fc G {I,..., 4, a G Ed} 

( 11 ) 

Here, E^{X) := Ek{PaX), where cr denotes an element of 
the permutation group Ed of d elements, and Per denotes 
the corresponding d x d permutation matrix, i.e. Pa-X = 
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(Acr(i), Acr(2), ■ ■ •, K(d)) = Act. TMs CRii bc easily seen by 
noting that for any X G R‘^ we have Ek{Pa\) < Ek{X^) 
for any k and any permutation ct € Sd- Hence, the ma- 
jorization conditions given in Eq. (4) are fulfilled iff the 
conditions in Eq. (11) are. 

For reasons that will become clear below, we consider 
the set of all Schmidt vectors that majorize A(^/’), not 
just the sorted ones, which is given by 

= {A s.t. A ^ A(^/>)}. (12) 

It is a well-known fact that (ip) is a convex poly¬ 
tope whose vertices are given by the vectors X{'ip)a-, with 
a € [20]. Hence, its H-representation reads 

= conv ({A(V')<T}aGEj ■ (13) 

Note further that, as the cardinality of is d!, every 
generic state in M.s{ip), i.e. every state with pairwise 
different Schmidt components, is represented d! times 
in This does not hold for non-generic states. 

However, as the set of non-generic states is of measure 
zero, we obtain that 

Vsipj) = ^i{Ms{ip)) = ^/r(-Adf'^(V'))- (14) 

We are going to show below that {ip) is generically 
a (d — l)-dimensional polytope. Hence, we choose for 
/i the Lebesgue measure in in order to arrive at 

a physically meaningful quantification of the source set. 
The non-generic case will be treated afterwards (see Ap¬ 
pendix A). 

In order to determine now the volume oi we use 

the results presented in [21]. There, the volume of a 
simple polytope, P C R^, with dim(P) = k, was shown 
to be 

m( 4^) = |det(m„)| ^ -■ (15) 

vevert{p) n(ei(^'),0 

Here, is the k x k matrix [ei{v)] e2(u);...; efc(r’)], 
where Cipv) = v — Vi{v) denotes the edge vector that 
connects the vertex v with its i-th neighboring vertex 
Vi{v) [30]. Moreover, the vector ^ is an arbitrary 

vector with the property that it is not orthogonal to any 
edge vector, i.e. {i,ei{v)) ^ 0 for all v and i. In fact, 
ldet(m„)j /k\ is the volume of the simplex defined by the 
vertex v and its neighbors. With the term simplex we 
refer to a fc-dimensional polytope in R^ that is the convex 
hull of fc -I- 1 vertices. The volume of the whole polytope 
P C R^ is thus the weighted sum of these different, but 
overlapping, simplices. 


Hence, in order to compute Vs{ip) via Eq. (15) it 
remains to show that {ip) is a simple polytope 

and to find the vector ^ with the properties mentioned 
above. Note that, due to the normalization condition, 
Ai = 1, we have that M^^{ip) is at most a (d— 1)- 
dimensional polytope in R'^. As a consequence, we need 
to identify the (d — l)-dimensional subspace of A, 
which contains A 4 g^{ip) and show that {ip) is sim¬ 
ple. More specifically, we are going to show that every 
vertex of Aipp^{ip) has generically (d — 1) neighbors and 
that the vertex and its neighbors are affinely indepen¬ 
dent, i.e. the (d — 1) edge vectors {ei{v)}i are linearly 
independent. This implies that Aipp^{ip) is a (d — 1)- 
dimensional simple polytope. Finally, we have to find ^ 
fulfilling the properties mentioned above. 

Let us first describe the (d — l)-dimensional subspace 
of RP^, which contains Alpp^{ip). Note that Alpp^{ip) C 
{A G R!^ s.t. Ylii Ai = 1} =: Al, where A = ^^ 

is a (d— l)-dimensional affine subspace of [31]. Here, 
!p+ = ^(1,1,...,!) and 

H = {A G 1?^ s.t. ^Xi = 0}. (16) 

i 

We consider hrst the generic case. That is, there 
exists no pair of different indices, i,j, such that Ai = Xj. 
Notice then that in this case the vertices of this polytope 
are in one-to-one correspondence with those of the so- 
called permutahedron, i.e. the polytope whose vertices 
are formed by all permutations of the coordinates of 
the vector (1,2, ...,d), which has been well-studied in 
polytope theory. It is known that the permutahedron 
is simple and its vertices have been characterized [19]. 
For the sake of readability, we give here a self-contained 
proof that Ai^^{ip) corresponds generically to a simple 
polytope in the (d — l)-dimensional subspace A, i.e. 
that every vertex has exactly d — 1 affine independent 
neighbors. 

First of all, we determine the neighbors of 
the vertex X{ip) = X'^{ip), which we assume 

to be sorted already. To do so, we consider 
here the set of d — 1 independent inequalities 
= {Ffe(A') < Ek{X{ip)) s.t. k G {l,...,d- 1}} 
and the set 122 = 12\12i, where 12 is as defined in Eq. 
(11). Clearly, X{ip) trivially satisfies all d — 1 inequalities 
in I2i as equalities. Its neighboring vertices must satisfy 
d — 2 inequalities in I2i and one independent inequality 
of 122 as equalities, as they must be elements of d — 2 
facets of which X{ip) is an element too. Using now 
that each of the vertices must be of the form X{ip)cr 
(see Eq. (13)), we have that the neighbors of X{ip) can 
only be of the form Ai,i = Hr-A('!/'), where 
denotes the permutation matrix permuting component i 
with component i + 1 and leaving the rest unchanged, 
e.g. Pti,2H'4’) = (A 2 , Ai, As,..., Ad). As there are 
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d — 1 possible permutations of this kind, we have 
d — 1 neighbors of the vertex [32]. As the d — 1 
edge vectors = X{tp) — X{ip)i are of the form 

e{X{ip))i = (0,..., 0, Ai+i - Xi,Xi- Aj+i, 0,..., 0), they 
span a (d — 1) dimensional space. Hence the polytope 
is indeed (d — l)-dimensional. The same 
argument holds for any other vertex A('0 )ct, which has 
the neighbors A^y = PaTVi.j+iA('0) [33]. Hence, all 
vertices have exactly d — 1 neighbors, which implies 
that the (d — l)-dimensional polytope is simple. Note 
that it does not need to be simple in case of degeneracy, 
i.e. for non-generic states. However, using continuity 
arguments, one can show that the volume of the cor¬ 
responding polytope can be computed as in case of no 
degeneracy (see Appendix A for details). 


We can now use the volume formula given in Eq. (15) 
in order to calculate the source volume of a generic 
state with Schmidt vector A(^/’) = X^{tjj)- The volume 
is translationally invariant, implying that = 

{tjj)) = ~ i/'/dff)^). The vertices of 

the translated polytope are given by X{'ijj)cr — 1/y/d^'^. 
Note that these vertices are d-dimensional vectors. How¬ 
ever, they obviously all belong to the (d— l)-dimensional 
subspace U defined in Eq. (16). In order to use Eq. (15) 
all vectors needed there must be considered as linear com¬ 
binations of an orthonormal basis of U, e.g. the matrix 
niy contains the coordinates of the edge vectors {ei{v)}i 
with respect to that basis. However, we use the follow¬ 
ing properties to circumvent this basis change. Note that 
the term jdet(TO„)j in Eq. (15) gives the volume of the 
/c-dimensional parallelotope that is defined by the edge 
vectors {ei{v)}i. Its numerical value is identical to the 
volume of a (fc -b l)-dimensional parallelotope of which 
the k vectors that define the base coincide with the orig¬ 
inal ones and the (fc -I- l)-th vector is normalized and 
orthogonal to the subspace that contains the base, e.g. 
the area of a rectangle is equal to the volume of a cuboid 
that is of height one and whose base coincides with the 
rectangle [34]. Eor the current computation this means 
that 


|det(TOAj| 


det 


-I-. 


; ei(Acr); e2(Acr);...; Cd-iiXc 


where as defined before, is the vector that is or¬ 
thogonal to the edge vectors ei{Xa) = Pae{X{ilj))i, with 
e(A('0))i = (0,..., 0, Ai+i - Xi,Xi- Aj+i, 0,..., 0). That 
is, we do not need to express the edge vectors in terms of 
an orthonormal basis of U in order to compute this de¬ 
terminant, but we can continue with the representation 
in terms of their d coordinates in the standard basis of 
R‘^. Using the fact that Pcr^'^ = and jdet(Pcr)| = 1 it 


follows that 
|det(mA,,)| = 

= ldet(P<,)l |det (<()+; ei(A('!/')); e2(A('0)); ■ ■ •; ed-i(A('!/'))) 

d-1 

fc=i 


The last expression can be easily obtained as each vector 
ei{X{ijj)) contains only two non-vanishing entries. It is 
now easy to choose a vector ^ G U fulfilling that it is non- 
orthogonal to any of the edge vectors ei(v), as required 
by Eq. (15). Note that ^ only appears in inner products. 
Moreover, the inner product is invariant under a basis 
change from an orthonormal basis of U, extended by </>+ 
in order to be a basis of i?'^, and the standard basis of 

. Consequently, we can again represent ^ and all other 
vectors in these inner products in the standard basis. In 
fact, for 5 = (1, 2,..., d) — ^^(I, !,...,!) £ U it holds 
that 

(e,e.(A.)) = {^,P.ei{X{m = (P<,-i^,e.(A(V>))) 

= (cr“\i) - a~^{i + 1)) (Ai+i - Xi) , (18) 

where we used that Pj = P^-i- Furthermore, we have 

(C, A(V')<t - 1/Vd0+) = Y1=i cr~\k)Xk - Plugging 
these expressions into the volume formula in Eq. (15), 
relabeling with ct, using that I4('!/') = ^{ip)) 
and the fact that A('0) is a sorted Schmidt vector, we 
obtain 


Vs{ip) 


1 ^fd 

d\yrY)\^ 

(Tt 


{ E cr(fc)Afc 

\fc=l _ 

d-1 



n cr(fc) — a{k -b 1) 

fc=i 


We show in Appendix A that this formula holds also 
in the case of degenerate Schmidt coefficients. Note 
that the separable state, \ip)g^p^ with Schmidt vector 
(1, 0,0,..., 0) can be obtained from any other quantum 
state via LOCC and therefore maximizes the source vol¬ 
ume, i.e. sup^g^ I4(()i) = Vs{ipsep)- According to Eq. 
(13) the vertices of the source set of |'0)sep are given 
by the standard vectors in R'^. {ipsep) is thus 

a simplex defined by the d vertices {ei}(Li whose vol¬ 
ume is easily computed. It is then straightforward to 
see that Vg{ipgep) = Vd/{d\{d — 1)!). Furthermore, it is 
easily seen that 14(^“*") = 0, as we show in Appendix 
A. Using the definition of the source entanglement, i.e. 
Es{ip) = 1 — we obtain the following lemma. 

Lemma 1. The source entanglement of a bipartite state, 
\ijj) £ QJ‘^ with sorted Schmidt vector X{ijj) is given 
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by 

/'dr ^ ^^ 

( E cr(fc)Afe - ^ ] 

EsW = 1 - E - - -• (20) 

o-6Sd _|_ 

fc=l 

In what follows we show how one can generalize the 
formula for the source entanglement in Lemma 1 in order 
to obtain a whole class of new operational entanglement 
measures as mentioned in Section II. 

Generalizations of the source entanglement 

We consider, for given \tjj) with sorted Schmidt 
vector X{ip), its source set of states |'I') € (F^ ®Q!^, 
where k > d. More specifically, we measure the set of 
states with greater or equal local dimensions that can 
be converted to |'0) via LOCC. This set is denoted by 
M^(^) and reads 

= {l^-) 00’'= s.t. I^-) IV')}. (21) 

In fact, we can identify |V’) with a state |d>^('0)) € 
(7^ 0 (7^ that has the Schmidt vector A(^'^(V')) = 
(Ai,..., Ad, 0,..., 0) G 1?^, where we simply appended 
k — d zeros to the initial, d-dimensional Schmidt vector 
of IV’). The volume corresponding to is then given 

by Vs if’) = 14(^^(V'))- In Ibis way we obtain a whole 
class of operational entanglement measures that are a 
generalization of the source entanglement presented in 
Lemma 1 and that read 

^'(^) = - k > d. (22) 

sup 

I0)e« 

Here, we divided by the supremum over all \(j)) GQ!^ 
such that the range of E^{'ip) is [0,1]. Note that these 
generalizations include the initial source entanglement 
for k = d. E^{iIj) is a polynomial of degree at most k—1 
in d — I Schmidt coefficients of IV”), where the remaining 
coefficient is given by normalization. We thus have ex¬ 
plicit formulae for a whole class of operational entangle¬ 
ment measures that quantify the set of LU-inequivalent, 
bipartite pure states with the same or higher dimension 
that can be transformed to a given bipartite state, |V’), 
via LOCC. Stated differently, the set of measures 

{E^mk>d. (23) 

characterizes how easy it is to generate a single copy of 
IV') from a single copy of another bipartite pure state via 
LOCC. 

C. Accessible volume 

We consider now the accessible set given in Eq. (8) 
for a state |V’) whose sorted Schmidt vector is given by 


A(V’) = Al(V’). As seen from Eq. (8) this set is the 
intersection of the halfspaces defined by the inequalities 

Ek{X')>Ek{\m Vfce{I,...,d-I}, (24) 

a; > A '2 > ... > A'rf > 0. (25) 

Our aim is now to compute the volume of this con¬ 
vex polytope. Let us first note that we consider here 
only the LU-equivalence classes as we are fixing the or¬ 
der of the Schmidt coefficients. Any other fixed order 
would obviously also represent all LU classes, despite the 
fact that the corresponding vectors belong to a different 
subset of Moreover, the volume of these different 
sets is the same. In contrast to the source set, they are, 
however, not connected in general and the union of these 
sets of differently sorted vectors that majorize A(V') would 
not lead to a convex set. In order to see that consider 
A(V’) and A' = (1,0,..., 0), i.e. the Schmidt vector cor¬ 
responding to a separable state. Clearly, it holds that 
X' >- A(V’). However, the set is in general not convex, as 
1/2A'-|-1/2 Pt-j 2 A' A if Ai > 1/2. The above mentioned 
properties of the union of the different representations of 
the accessible set is the reason why we do not consider 
this set, although we used such an approach in case of 
the source set. However, the accessible set defined by the 
halfspaces in Eq. (24) is, as the source set was, a convex 
set. 

Let us now demonstrate how the volume of A4a(V’) can 
be computed. To do so, we first determine the vertices 
of the accessible set. Thereby, it is useful to get rid of 
the last Schmidt coefficient, which is determined by the 
normalization condition, i.e. A^ = 1 — Ej^_i{X). Hence, 
from now on, we consider the Schmidt vectors as d — 
1 dimensional vectors and the ^-representation of the 
accessible set is given by 

Afc(A') > Efe(A(V’)) Vfce {!,...,d-1}, (26) 

A;>...>AEi>l-A,i_i(A')>0. (27) 

A vertex of Ma, Vi, is the unique vector in Ma C , 
where d— 1 independent inequalities are satisfied as equa¬ 
tions. Let us remark here that the method to calculate 
the vertices of the accessible set and its volume differs 
from the one used to calculate the source volume, as we 
will explain in the following. To determine the source 
volume, we calculated the (d — I)-dimensional volume 
of the source set that is contained in A G (respec¬ 
tively U, after a translation). This is why we had to, 
at least conceptually, transform into an orthonormal ba¬ 
sis of U. To calculate the accessible volume, however, 
we abandoned the last Schmidt coefficient. As a result, 
the id-representation in Eqs. (26) and (27) describes the 
projection of the accessible volume that is contained in 
A G R^ onto the subspace spanned by the first (d — 1) 
standard basis vectors. Subsequently we identified this 
subspace with R^~^. As shown in Appendix B, the re¬ 
sulting volume is 1/i/d-times as big as the original one. 





However, this constant factor is at the end irrelevant as 
we rescale the volume in order to obtain the accessible 
entanglement (see Eq. (1)). 

Let us now consider the determination of the vertices 
of the accessible set. Note that X^ip) is always a vertex as 
it obeys all d — 1 inequalities in Eq. (26) as equalities. As 
we have 2d— 1 inequalities, there can exist at most 
vertices. Clearly, there will be much less in general. Con¬ 
sider for instance for d = 4 the state 1/100- (30, 27, 24,19) 
that has ten vertices, while the state 1/10 -(4, 3, 2,1) has 
only eight, which is a lot less than (g) = 35 . As can be 
seen from these examples, even the number of vertices 
depends also for non-degenerate states very strongly on 
A(^). Furthermore, it is not clear if the accessible set 
is always a simple polytope for non-degenerate states. 
These are probably the reasons why it is difficult to de¬ 
rive a closed expression for the vertices or the volume of 
the accessible set given some state corresponding to the 
Schmidt vector A(^/’). The vertices can however be easily 
computed for a given X{^p), as we will briefly explain in 
the following. 

Assume that one vertex, vi, is known (e.g. X{ip) in the 
case of the accessible set). Let denote a set of d — 1 
independent inequalities which are fulfilled as equalities 
for the given vertex. In order to find a neighboring vertex, 
pick one inequality which is in l7i and replace it by an 
inequality from the id-representation which is not in I7i. 
The vector in R‘^~^ satisfying these inequalities now as 
equalities (and of course obeying all other inequalities) is 
a neighbor of vi , as it is an element of d — 2 facets that 
also contain vi. One can continue in this way until one 
finds all possible vertices. 

In the literature there exist several algorithms for com¬ 
puting the vertices and the volume of arbitrary convex 
sets. The reader is referred to [22] for a review on the 
basic properties of some algorithms for volume compu¬ 
tation. Usually the volume is computed by using either 
the so-called triangulation method or the signed decom¬ 
position method. In the former method the polytope is 
decomposed into simplices with mutually disjoint interior 
such that the volume of the whole polytope is the sum 
of the volumes of the individual simplices. In the latter 
method the simplex is decomposed into signed simplices, 
i.e. simplices to which a positive or negative sign is asso¬ 
ciated, that are allowed to overlap. These are then, de¬ 
pending on their sign, added or subtracted successively 
in order to obtain the volume of the polytope. We will 
use here the algorithm presented in [23] to calculate all 
vertices of the accessible polytope. Although one can find 
the vertices using the method outlined above, the aim is 
of course to perform the computation using few resources. 
In [23] it has been shown that, by choosing the inequali¬ 
ties that have to be interchanged as explained above ac¬ 
cording to certain rules, the vertices can be found easily. 
In fact the algorithm finds the N vertices of a polytope in 
R‘^~^ defined by a nondegenerate system of 2 d —1 inequal¬ 


ities in time 0{d^N) and 0{d?) space. Note, however, 
that N can in principle grow exponentially with d. A 
revised version of the algorithm presented in [24] is also 
capable of computing the volume of the convex polytope 
using the triangulation method described above. 

Whereas the numerical methods described above can 
be utilized to compute the accessible volume for arbitrary 
states, we present in Sec. Ill D some examples of bipar¬ 
tite systems of small dimension for which we computed 
analytic expressions for the vertices, the corresponding 
accessible volume and the associated accessible entangle¬ 
ment. Note that, similarly to the source entanglement, 
also the accessible entanglement can be generalized by 
considering different Hilbert spaces of the initial state 
and the state in the accessible set ([14], see Sec. HID for 
examples). 

Before concluding this subsection, let us present some 
states which always correspond to vertices of the ac¬ 
cessible polytope. The reason for doing so is that the 
knowledge of the vertices of the accessible set allows us 
to gain more insight into the entanglement properties of 
the state. The following vectors are always vertices of 
the accessible set of a state with sorted Schmidt vector 
A(V-) = (Ai,...,Ad), 

V\ — (Ai,Ai,...,Ai, Xnorm,l ib;---;!!)! (^^) 

V2 = (Al, A 2 , A 2 , . . . , A 2 , Xnorm,2, 0 , . . . , 0 ), (29) 

^(1-2 = (Al, A 2 , A 3 , . . . , Ad- 2 , Ad- 2 , A nor771,cZ— 2 ) ; (30) 

where 0 < Xnorm,i < Xi is such that all components add 
up to one. The reason for that is that Vi obviously ma¬ 
jorizes A(V') and is an element of d— 1 facets of the acces¬ 
sible set. More precisely, Vi satisfies i of the inequalities 
from Eq. (26) and d — 1 — i inequalities from Eq. (27) 
as equalities. Note that these d — \ inequalities are inde¬ 
pendent such that Vi is indeed a vertex. Furthermore, all 
maximally entangled states of two A:-level systems, ) , 
i.e. states where the first k entries of the sorted Schmidt 
vector are 1/k and the remaining ones are zero, are in 
the accessible set iff Ai < 1/k. To see this, note that this 
criterion is equivalent to Ei{(j)//) > Ei{X{tj})). Since A('0) 
is sorted, this implies that Ei{(j)//) > Ei{X{'ip)) for all i 
and therefore (see Eq. (4)) that is accessible from 
ji/). It is straightforward to show that X{(j)//) is in this 
case also a vertex of the accessible polytope that fulfills 
d — 1 of the inequalities in Eq. (27) as equalities. 

D. Source and accessible volume for low 
dimensional bipartite systems 

In this section we use the results obtained above in 
order to compute the accessible as well as the source en¬ 
tanglement of bipartite systems. More specifically, we 
present the source and accessible set of small dimensional 


9 


systems, i.e. d=2,3,4inFigs. 2-4 and briefly discuss 
their properties. By means of these figures we can also 
derive closed expressions for the accessible volume of a 
state of two qubits or two qutrits. 


1. Two qubits 

In the simplest example of a two-qubit state there ex¬ 
ists only one independent Schmidt coefficient, such that 
the source and accessible volume of a state with sorted 
Schmidt vector A(^) = (A, 1 — A) read 

MsW = {(A'l, 1 - A'l) s.t. 1/2 < a; < A}, 

MaW = {(A'l, 1 - A'l) s.t. A < a; < 1}. 

Consequently, the two sets are complementary in the 
sense that their union results in the whole state space 
and there do not exist LOCC-incomparable states. That 
is, given two states one of them can always be converted 
into the other via LOCC. For this simple example the 
volumes can be directly read off from Figure 2, in which 
both sets are depicted for a given state. It is easy to 
see that the source and accessible volume of a two-qubit 
state lip) are given by 

Vs{iIj) = V2{X,-1/2), 

VaW = V2{l-Xi). 


Using these formulas it is straightforward to see that the 
source and accessible entanglement (see Eqs. (1)) of a 
two-qubit state are in fact identical and given by 

Esi^) = Ea{iP) = 2(1-Xi). 

Clearly, each of them gives a unique characterization 
of the entanglement contained in a two-qubit state, i.e. 
they allow to recover its Schmidt coefficients. 


2. Two qutrits 

Let us proceed with the investigation of a system of two 
qutrits. In contrast to the two-qubit case one can find 
two-qutrit states that are LOCC-incomparable. Conse¬ 
quently, the source and accessible set of a generic two- 
qutrit state are not complementary, as can be seen in Fig. 
3. More precisely, their union does not lead to the whole 
state space. Moreover one can easily see that represen¬ 
tations of the accessible set that correspond to different 
orderings of the Schmidt coefficients are not connected 
to each other. In contrast to that the union of the corre¬ 
sponding representations of the source set, i.e. A4f'^('i/'), 
is a convex set. We use Eq. (19) to obtain the source vol¬ 
ume of a state \tjj). Based on Fig. 3 we can also calculate 


A 2 



FIG. 2: (colors online) The source and accessible set of a two- 
qubit state with sorted Schmidt vector A = (0.6, 0.4). The 
line that connects the points ei = ( 1 , 0 ) and 62 = ( 0 , 1 ) cor¬ 
responds to the set of valid, i.e. normalized, Schmidt vectors. 
The shaded region highlights the set of sorted Schmidt vectors 
that is in one-to-one correspondence to the LU-equivalence 
classes. The red (blue) lines in this region depict the source 
(accessible) volume of the quantum state respectively. 


its accessible volume using elementary geometry. That 
is, we subdivide the accessible set into mutually disjoint 
triangles. The shape of the accessible set, and with it the 
triangulation, depends on the largest Schmidt vector, Ai, 
of the state. This can be easily seen in Fig. 3. The vol¬ 
ume of a triangle with vertices vi , V 2 and V 3 is then given 
by the formula 1/2 |(?;2 — 'Ci) x (^3 — ?;i)|, where a x b 
denotes the cross product between vectors a,b £ R^. In 
analogy to the generalization of the source volume dis¬ 
cussed in Sec. IIIB we moreover compute I//('0), i.e. the 
volume of accessible states with at least one zero Schmidt 
coefficient. Stated differently, we compute the volume of 
two-qubit states that is accessible from lip). We obtain 
the volumes 


VsW 

Vaitp) 


6 A 2 A 3 — 3 A 2 -l- 6 (A 3 — 1)A3 -I- 1 

4v^ 



A 2 A 3 — - (1 — 2Ai)^ 


ifAi>i 

ifAi<i 



Ai) 


if Ai < f 


Starting from a state \ip) with three nonzero Schmidt 
coefficients one can, as pointed out in Sec. IIIC, access 
the maximally entangled state of two-qubits, |(/^), with 
Schmidt vector A(^ 2 ) = (1/2,1/2,0), iff Ai < 1/2. It 
is clear that in this case one can obtain any two-qubit 
state. This is why such a state maximizes V^{'ip). It is 
straightforward to compute the entanglement measures 
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from the volumes above. We obtain 


3. Two four-level systems 


— 3A2 — 6A2A3 — 6(A3 — 1)A3, 


Ea{tp) 

ElW 


( 12 A 2 A 3 

[i2[A2A3-1/4(1-2Ai) 2] 
f 2(1 — Ai) if Ai > - 


ifAr>i 

ifAi<i 


where we defined EKip) = sup^v^\,p) - Using Eq. (22) 
one can furthermore calculate 


Finally, we also consider quantum states of two four- 
level systems. Also for such a system we can depict the 
source and accessible set of a given state (see Fig. 4), 
as they are at most three-dimensional. We could not 
give an analytic formula for the accessible entanglement. 
However, the algorithm described in [24] (see also Section 
III C) can be used to find the vertices and the volume of 
the accessible set of a given two-qutrit state. We used 
this algorithm to determine these properties for the state 
depicted in Fig. 4. Analogously to before, the source 
entanglement can be computed using Lemma 1. It is 
given by 


EtW = {2X1 + BA^Aa + 3(3 - 4A2)A2 - lOA^) . 

Note that one can show that Egiip) together with E'^i^ji) 
uniquely characterize the entanglement contained in a 
quantum state of two qutrits. That is, these two mea¬ 
sures uniquely define the Schmidt components of the 
state. Hence, its entanglement is completely character¬ 
ized by how many two-qutrit states and states of two 
four-level systems can reach the state at hand. 



FIG. 3: (colors online) (a) The source and accessible set 
of a two-qutrit state with sorted Schmidt vector A = 
(0.6,0.37,0.13). The triangle with the three vertices ei = 
( 1 , 0 , 0 ), 62 = (0,1,0) and 63 = (0,0,1) corresponds to all 
valid, i.e. normalized, Schmidt vectors. The thick, dashed 
line encloses the set of sorted Schmidt vectors that is in one- 
to-one correspondence to the LU-equivalence classes. The red 
(blue) regions depict the source (accessible) set of the quan¬ 
tum state respectively. The hyperplanes that correspond to 
the halfspaces defining the set {ip) are also indicated by 
the dotted lines. The white dot corresponds to the maxi¬ 
mally entangled state \(p 2 ) of two qubits with Schmidt vector 
(1/2,1/2, 0). Since the depicted state has Ai > 1/2 the state 
\<Pt) 1® element of its accessible set. Its accessible set 

has four vertices. 

(b) The source and accessible volume of a two-qutrit state 
with sorted Schmidt vector A = (0.47, 0.36, 0.17). In contrast 
to the state in (a), this state fulfills Ai < 1/2 s.t. the source 
volume contains the state j (p 2 ) and the accessible volume has 
five vertices. 


Es{ip) =4A^ -h 12A^A3 - 24A^A4 - 24A2A^ -h 24 A 2 A 3 A 4 
-h 12A2A^ - 20Ai -h 12A^A4 -h 18A^ -h dSAgA^ 
- 36A3A4 -h 20A^ - 3OA4 -h I2A4. 


A 2 



FIG. 4: (colors online) The source and accessible set of a 
state of two four-level systems with sorted Schmidt vector 
A = (0.4, 0.3, 0.2, 0.1). The valid, i.e. normalized, Schmidt 
vectors correspond to all states in the set {(Ai, A 2 , A 3 ) s.t. Ai > 
X)i=i Ai < 1}. The last coefficient is given by normaliza¬ 
tion. The convex set indicated by the shaded region (in¬ 
cluding the red and the blue region) highlights the set of 
sorted Schmidt vectors that is in one-to-one correspondence 
to the LU-equivalence classes. The red (blue) regions depict 
the source (accessible) set of the quantum state respectively. 
The vertices of the accessible set and its volume have been 
computed with the algorithm presented in [24] (see also Sec. 
IIIG). Weobtain Ea{ip) = 0.904 and Ea{ip) = 87/125 « 0.696 
for the depicted state. 


IV. ENTANGLEMENT OF MULTIPARTITE 
SYSTEMS 

In the previous section we showed that one can obtain 
explicit formulae for a whole class of new operational 
entanglement measures that characterize how easy it is 
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to generate a single copy of a quantum state from a single 
copy of another quantum state via LOCC. Moreover, in 
[14] we proved that these measures can also be used to 
completely characterize three-qubit entanglement. Here, 
we will determine the source and accessible entanglement 
of generic four-qubit states. 

Clearly the source and accessible volume can only be 
computed if one knows all the possible LOCC transfor¬ 
mations of the states of interest. Note that in contrast 
to the bipartite case, the characterization of all possi¬ 
ble LOCC transformations is generally very difficult, as 
also protocols involving infinitely many rounds of com¬ 
munication have to be considered [15]. However, in this 
section we will derive the necessary and sufficient condi¬ 
tions for the existence of a LOCC transformation among 
pure states by completing the results on deterministic 
state transformations of four-qubit states presented in 
[9]. This criterion allows us then to compute the two en¬ 
tanglement measures. Note that we want to consider de¬ 
terministic transformations among fully entangled states, 
which are only possible among states within the same 
SLOCC class. In [25] it was shown that there exist in¬ 
finitely many SLOCC classes for four-qubit states. Here, 
we only consider generic four-qubit states, which belong 
to the SLOCC classes denoted by Gabcd with representa¬ 
tives of the form [25] 

\'^)seed = ^ (| 0000 )+| 1111 )) + ^ (| 0011 )+| 1100 )) 

+ ^ (| 0101 )+| 1010 )) + ^ (| 0110 )+| 1001 )), 

(31) 

with b,c,d G(F, a & R and + \b\^ + jcj^ -I- \df‘ = 1, 
6 ^, c^, b^ d^ b^ and the parameters ful¬ 

fill the condition that there exists no q gQJ\{1 } such that 
{ } = { qa^ , qb"^ , qc^ , qd^ }. These states are re¬ 

ferred to as the seed states, which are parametrized by 
six real parameters, due to normalization and the global 
phase. Note that Py <^1 \'^)seed for any 

particle permutation Py with i ^ j, i,j € { 1 , 2 ,3,4}. 
That is, any seed state has the property that permut¬ 
ing the particles leads to a LU-equivalent state. Any 
state in some SLOCC class Gabcd can be written as 
|4'((;)) = g \^)gaed^ a local invertible operator (not 
necessarily determinant 1 ) g G I/, i.e. g = 

G GL{2), and the corresponding seed state \^) gf.ed- 
Ignoring for the moment the normalization, the positive 
operators G* = (g*)fg* fulfill without loss of generality 
tr(G®) = 1. Furthermore, they are positive full-rank op¬ 
erators written as G* = 1/21 -I- with 7 ^ G P 

and 0 < [ 7*1 < 1/2, where 7 * = ( 71 , 72 , 73 )- Sometimes 
we will also use the notation 7 [„ with w G {x, y, z} for 
the state parameters. Note that the symmetries of the 
seed states are given by {crf'^jLo. i-e. af'^ \'^)seed = 
Iseed’ *• These symmetries can only simulta¬ 

neously change the sign of two parameters 7 ^ and 7 ^, 


with k,l G {1,2,3}, k ^ I, for all i. Thus, G* can be 
made unique [9]. By choosing = VG^ and sorting the 
parameters of the seed states the form 

|^(fl^5^d^5^)) =fll^).eed (32) 

is a unique standard form of the state {g^, g"^, g^, g‘^)y 
Therefore, generic four-qubit states are in the same LU- 
equivalence class iff their standard forms coincide. As we 
only consider representatives of LU-equivalence classes 
for the computation of the source and accessible volume 
we will pick the state \'^{g^,g^, g^, g^)) as such. 

In [9] the maximally entangled set (MES) for the generic 
four- qubit states has been derived. The MES is the min¬ 
imal set of states from which any other genuinely mul¬ 
tipartite entangled state can be obtained via LOCC. It 
has been shown that this set is of full measure for four- 
qubit states, in contrast to two- and three-qubit states. 
The reason for that is that most four-qubit states are 
isolated, meaning that they can neither be reached nor 
converted into any other non LU-equivalent state. Hence, 
for isolated states the source and accessible volume van¬ 
ish. The non-isolated states have been shown to be (up 
to permutations of the local operators) [9] 

= 9^ ®gl®gl®9t,\^)seed^ (33) 

with gl, G span{l, cr^}, such that G[„ = {gl^Ygl, = 
1/21 -I- ^IjCTw, w G {x,y,z} and arbitrary g^ as defined 
above. Note that we will in the following choose the two 
parameters 7} and 7} nonnegative and allow for negative 
values of the third parameter 7}, i.e. 7}, 7} G [0,1/2), 
7{ G (—1/2,1/2). We can define the standard form of 
\'^{9^,9i,9l,9t)) like that, as the symmetries of the 
seed state can always change the sign of two parame¬ 
ters of 7 ^ simultaneously as mentioned above. Hence, 
the non-isolated states form a 12 -parameter family. As 
only this zero measure set of states allows for non-trivial 
LOCC transformations these are the only states we need 
to consider here. Note that we consider here states up to 
permutations of the parties as they do not alter the nec¬ 
essary and sufficient conditions for LOCC convertibility 
and hence, we can easily compute the source and acces¬ 
sible volume of the permuted states (see below). The 
non-isolated states in the MES, that can not be reached 
by any other state but can access other states via LOCC, 
form a 10 -parameter family with g = glu ® at) ® g% ® gtj^ 
w G {a:, y, z} (excluding 9 ^ 1 for exactly one i) [9]. 

In the following we will derive the necessary and 
sufficient conditions for LOCC transformations among 
generic four-qubit states. Subsequently, we will always 
denote the state of interest by \'^{g^, 9 ^, 9 ^, g'^)) as in 
Eq. (32) and the states in the source or accessible set 
of the state of interest by |$(/i^, , h^)) = h ['1')^,^^^^. 

Note that one obviously only has to consider non-isolated 
states. Eirst, one can easily show that LOCC transfor¬ 
mations exist only among very particular pairs of states, 
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as stated in the following observation which is proven in 
Appendix C. 

Observation 2 . A state |'I'(g)) can be transformed to 
a reachable state h^)) by LOCC only if 

9 = 9^ ® 9w^ 9w^ 9w! same w € {x, y, z}. 

Note that not any state of the form h^, 

is reachable, as stated in the next lemma, which has been 
proven in [9]. 

Lemma 3. A generic non-isolated four-qubit state 
|<l)(/i^, h^, h^, /i^)) = h can be reached via LOCC 

by some other state \'i', 9w t 9 w^ 9w): 
iff (up to permutations) either 

(i) h = h^ ® hfj ® h)^ ® h)j, with h^ ^ and at 

least one i € {2,3,4} such that pf 1 and h 

h} 0 0 1 0 1 for V w, V £ {a;, y, z} or 

(a) h = ft,} 0 (g) 1 (g) 1, with ft}, pt 1, ft} pt 1 and 

V w, V € {x, y, z}, or 

(Hi) ft = ft^ 010101 , with h^ pf 1 . 

As before we call |$(ft^, ft},, ft},, ft},)) with ft defined as 
in one of the three above cases a reachable state. Note 
that the reachable states are of course the non-isolated 
states excluding the non-isolated states in the MES. 

Let us now show from which states these states are 
reachable. In order to do so, let us denote by 'y\ the 
real state parameters of g},,g},,g},)), e.g. = 

= 1/21 + by the state parame¬ 

ters of |4>(ft^, ft},, ft},, ft},)) (with corresponding parame¬ 
ter vectors 7 * and Q respectively). Furthermore, we will 
use the notation 7 }(,„) (and the same for }}(„)) whenever 
we consider both the u and v component of 7 ^ for u u, 
{1,2,3}. 

Lemma 4. A state , 9^, 9^, gti)) can non-trivially 
access a reachable state |$(ft^, ft},, ft},, ft},)) within case 
(/) of Lemma 3 iff the following condition (for reach¬ 
able states) holds for L € {ftii,*^}. 

(i) C}, 7 ^ 0 for 

(a) at least two i € {1,2,3,4}, 7 }, = for 
3 = 1 , 2 ,3,4, and 7 }(„) = ( 2 p - 1 )C}(„) with 
p € [1/2,1) and {v,w,u} = {1, 2, 3}. 

(b) exactly one i G {2,3,4} and }}, = 0, C}(„) 7 ^ 0 
for {u,w,u} = {1,2,3} and either the condi¬ 
tions from above are fulfilled or 7 }, = 0 for all 
kf^i,k€ { 1 , 2 ,3,4}, 7 }(„) = 0 and 7 }, < Cw 

(a) 73 = 74 = 0 , 7 },(,^) = 0 , 7 }, < C and 7 } < C}, 
with {r, w, u} = { 1 , 2 ,3}. 

(Hi) 7 }, = 0 for i = 2,3,4 and either 


Cfc 7 ^ 0 G {1, 2,3}. Then with r^ = 7fc/Cfc ihe 
state parameters have to fulfill the three inequal¬ 
ities 1 -I- ri — r 2 — r 3 > 0 , 1 — ri -|- r 2 — ra > 0 
and 1 — ri — r 2 + r 3 > 0 , with 0 < r 3 ( 2 ) < 1 cind 
-1 < ra < 1. 

(b) Or Cl = 'll = 0 for exactly one k. Then Cl > ll 
and Cj > 7 } with {i,j, k} = { 1 , 2 , 3}. 

(c) Or cl 0 and thus, 7 } 7 ^ 0 for exactly one k. 
Then Cl > 7 } with k G {1,2,3}. 

The proof of Lemma 4 can be found in Appendix C. 
Note that Observation 2 together with Lemma 3 and 
Lemma 4 constitute the necessary and sufficient condi¬ 
tions for LOCC convertibility among generic four-qubit 
states. We will use them in the subsequent subsections to 
measure the entanglement contained in these states and 
summarize them now in the following table. We choose 
from now on without loss of generality w = x. 


Initial state 

Final state 

Necessary and suf¬ 
ficient conditions 

Case 

\ 4 ’{g^,gl,gl,gt)) 


7} = cl and 3 p G 
[1/2, 1) s.t. (2p - 

1 )C 2 ( 3 ) = 72 ( 3 ) 

(ia) 


|<f.(ft}(,),ft},l,l)) 

7i < Cl and 7^(3) < 
C 2 ( 3 ) 

(ii) 

1,1,1)) 

|cf>(fti, 1,1,1)) 

9 ^ 9 ^ l,ft}, ii,j,k} = 
{ 1 , 2 , 3 } and either 

• Ci 7^ 0 and 1 — 

Vi 3 - rj — Tk > 0 
with n = 7}/C}, 

• or =0 for 

exactly one k, 

7 fc = 0 , Ci > 3} 
and cj > yj, 

• or cl A 0 and 

7} / 0 for ex¬ 
actly one k and 

c} > -yl 

(iiia) 

(iiib) 

(iiic) 

| 4 '((/}, 1,1,1)) 

|cf>(fti, 1,1,1)) 
4 >(ft},ft},„l,l)) 

IV IV 

(iiia) 

(ib) 

l’I')seed 

4 >(ft\ 1,1,1)) 
4 >(ft},, ft},„,!,!)) 


(iiia) 

(ib) 


TABLE 1: Summary of all possible nontrivial LOCC transfor¬ 
mations of generic four-qubit states (up to permutations) with 
appropriate necessary and sufficient conditions (see Lemma 
4). Note that also the normalization condition on the pa¬ 
rameter vector, i.e. 0 < | 7 ^| , jC^j < 1/2 has to be fulfilled 
together with the necessary and sufficient conditions in the 
third column of the table. In the fourth column of the table 
the corresponding cases of Lemma 4 are given for each of the 
LOCC conversions. 

Due to the fact that the criterion for the existence of 
a LOCC protocol does not depend on the seed parame¬ 
ters, it is evident that the measures will only depend on 
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the state parameters and not on the seed parameters. To 
be more precise, we denote by |\E'(a, 6, c, the seed 

state given in Eq. (31). Then, the source and accessible 
entanglement of a state ® g'^ ® ® c, d)) 

coincides with the corresponding entanglement of the 


state g^ ® g"^ ® g^ ® g^ d>(a, b, c, d) 


for any choice 

seed 


of the seed parameters a, b, c, d and a, b, c, d. 


A. The source and accessible volume of generic 
four-qubit states 


(34) we get for an y reachable state that 
and 72 ^ < C 2 < a/(1/4 - + {llhlY) (72 can 

be chosen nonnegative in the standard form). Hence, the 
one-dimensional volume is determined by a line integral 
which results in 


Va{\'^{g^,gl,gl,gt)))=\Ji/^-{ilY-\J{nlY+{ilY- 

(35) 


Given the necessary and sufficient conditions for 
LOCC convertibility of the previous section, we will now 
compute the source and accessible volume of a state 
= 9\^)seed- Recall that we consider 
only fully entangled four-qubit states in the source and 
accessible set of Note furthermore 

that we will use the freedom of choosing different mea¬ 
sures g for computing the source and accessible volume, 
as explained in Sec. II. Whenever, the source or acces¬ 
sible states of certain states live in manifolds of different 
dimensionality, we choose different measures to compute 
their volumes (see also [14]). Otherwise, we would assign 
zero values to the volumes of some states, even though 
they can be reached or accessed by other states. Hence, 
by using different measures for the computation of the 
volumes we can compare the relative strength of states 
whose volumes have the same dimension and moreover, 
regard e.g. states with a two-dimensional accessible vol¬ 
ume obviously as infinitely more powerful than states 
with a one-dimensional accessible volume. We choose, 
as stated before, without loss of generality w = x, i.e. 

= 9^ ®9l®9l®9t\'^)seed^ as all the 
other cases can be treated analogously. We will in the 
following consider the computation of the volumes of 
states, for which different necessary and sufficient con¬ 
ditions have to hold, in the same order as in Table I. 
Hence, we start by computing the source and accessible 
volume of a generic state \^{g^,g^,g^,g'^)) with 7 ^ 7 ^ 0 
for at least two i G {1,2,3,4} (see Lemma 4 [case (ia)]), 
which is not in the MES, i.e. g^ 7 ^ (?{. As before we 
denote by 7 ^ the parameters of g^ and by the param¬ 
eters of . Due to Lemma 4 [case (ia)] the accessible set 
of parameters = (Ci) C21C3) is given by 

{C^ s.t. |C^| < 1/2 , Cj = 7 i\ (34) 

and 3p, 1/2 < p < 1 s.t. (2p - 1 )C 2 \ 3 ) = 72 ^ 3 )}- 

Note that the first condition on the norm of the vec¬ 
tor is due to the fact that has to be posi¬ 
tive semidefinite. Therefore, the accessible volume of 
with 7 { 7 ^ 0 for at least two i G 
{1,2, 3,4} is a one-dimensional line as the two param¬ 
eters C 2 s-nd can only be increased by the same pro¬ 
portion (see Fig. 5). Then from the conditions in Eq. 


Thus, with = 1/2 we get for the accessible entangle¬ 
ment Eai\'^{g^,gl,gl,g^))) = 2K(|^'(ff\ff^,5x>5x)))- 
Note that we obtain here and in the following the supre- 
mum of the accessible (source) volume by simply opti¬ 
mizing the corresponding volume, which is in this case 
given by Eq. (35), over the valid parameter space. Ob¬ 
viously, we get similar expressions for the accessible vol¬ 
ume of the states given by permutations of the parties 
of e.g. for the state \'i/{gl,g^,gl,gi)) 

the accessible volume is the same as in Eq. (35) with 7 / 
replaced by 7 ^ for fc = 1, 2, 3. 

Let us now compute the source volume of the state 
\'^{9\9l,9l,gt))- Due to Lemma 4 [case (ia)] 
the parameters of all states in the source set of 
with 7 { 74 0 for at least two i G 
{ 1 , 2 ,3,4} (the other cases are treated below) are ele¬ 
ments of the set 

{C^ s.t. 7 / = Cl and 3p, 1/2 < p < 1 s.t. 

(2p-1)72\3) = C2\3)}- (36) 

It is easy to see that |C^| < | 7 ^| and thus, |C^| < 1/2 
is fulfilled for all vectors Hence, the source volume 
of \'i'{g^,gxTgxT9x)) i® again a one-dimensional line de¬ 
fined by C 2 < 72 and the other parameter is fixed by 
C3 = C2 73/72- Therefore, we obtain the following source 
volume 


Vsi\'i’ig\gl,gl,gt))) = \J {ilY + inlY ■ (37) 


The corresponding source entanglement measure is with 
ysup = 1/2 given by Es{\'l!{g\glgl,gt))) = 1 - 
2Vs{\^{g^^g^^g^^gt)))- Note again that we get similar 
expressions for the source volume of the states that are 
given by permutations of the parties of |4'(p^, p^)) 

by simply replacing 7 / with 7 ^ for i G {2, 3,4}. 
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FIG. 5: (color online) Source (red line) and accessible (blue 
line) sets, Ms and Ma (see Eqs. (35) and (37) for the vol¬ 
umes), of a generic 4-qubit accessible state |'5((?^, 5^, 5x)) 

with parameters 71 = 0.15, 72 = 0.2 and 73 = 0.1. The 
shaded region corresponds to valid states in the standard form 
which are in one-to-one correspondence with LU-equivalence 
classes. 

Let us now consider states of the form 
\^{9x^9'xt9xi9x))^ where at least two state param¬ 
eters are non-vanishing, i.e. non-isolated states in the 
MES that are not the seed state. For these states the 
last conditions of Lemma 4 [case (ia)] are fulfilled by 
choosing p = 1/2, as 72 = 73 = 0. Hence, a state in the 
MES can reach any state of the form \^{h^ 
as long as 7/ = and |C^| < 1/2. Moreover, the 
local operator of each of the parties can be changed, 
that is also for instance the state |$(g^, 5;^)) 

can be reached as long as the corresponding conditions 
mentioned above are satisfied. Hence, the accessible set 
of parameters is given by 

{C s.t. Id < 1/2 and 7 i = Civ z G {1,2,3,4}}. (38) 

Therefore, the accessible volume of \'^{g].,g1,g^,g'^)) is 
the sum of four discs of radius Ri = ^1/4 — (7i)^i for 
i G {1,2, 3,4} i.e. 

Va{\'i>i9i,9l,9l,9t))) = 7r((l/4- (7i)^) + (1/4 - ( 7 ?)^)-^ 

(l/4-(7iT) + (l/4-(7i"d)- 

(39) 

It can be easily seen that the supremum of the acces¬ 
sible volume is given by = tt. Thus, we obtain 

for the accessible entanglement of convertible states in 
the MES Eai\^iglglgl,gt))) = (1/4- (yi)^) + (1/4- 
(7i)^) + (1/4— (7i)^) + (1/4— (7i)^)- The source volume 
of states in the MES is zero, as no state can reach a state 
in the MES via LOCC, i.e. Vsi\'^igl,gl,gl,gt))) = 0 
and hence, the source entanglement is equal to one, i.e. 
^^(|4'(5x,5x>5x.5x))) = 1- 

Let us now consider the source and accessible vol¬ 
ume of states as in Lemma 3 [case (ii)]. Here, we 
choose without loss of generality states of the form 
|4'(ffy,ff^,l,l)) Due to Lemma 

4 these states can only be transformed into states of 


the form |<i)(hj, 1,1)) (see also Table I, second row). 

The necessary and sufficient conditions for deterministic 
transformations into such a state are given in Lemma 4 
[case (ii)], which imply that the accessible set of param¬ 
eters is 


{(C2>C/) s-t. 72 < Cl and 71^ < C?} ■ (40) 

Hence, the accessible volume is equal to the area of a 
rectangle (see Figure 6), i.e. 

Vail'Siglgl 1,1))) = (1/2 - 7l)(l/2 - 7?). (41) 

Then, with 1 //“p = 1/4 the corresponding entan¬ 
glement measure is equal to Ea{\^{gyTg1^1,t)')) = 

414(|4'(gy,g^, 1,1))). For the source volume of 
|4>(d, 1,1)) the conditions on the parameters of 

states of the form |$(/iy, 1,1)) in the source set are 

given by 


{(Cl, Cf ) s.t. 7I > Cl and > C?} , (42) 

and thus, the volume is again equal to the area of a rect¬ 
angle (Figure 6), i.e. 

Vs{\'f{gl,gl,l,t))) = 4-fljf. (43) 

The source entanglement is given by 
Esi\'^igl,gl,tA)}) = 1 - Vsi\'^igl,gl,l,l))), as 

ysup ^ 2. 



FIG. 6: (color online) Source (red rectangle) and accessi¬ 
ble (blue rectangle) sets, Ms and Ma (see Eqs. (41) and 
(43) for the volumes), of a 4-qubit accessible state of the 
form [4'(y^,3^, 1, 1)) with parameters 71 = 0.1, 72 = 0.3. 
The dashed line encloses the grey area of valid states in 
the standard form corresponding to LU-equivalence classes. 
Hence, only states inside this grey area, where each of them 
is a representative of one LU-equivalence class, are in the 
source/accessible set. 


Next we examine the source and accessible volume of 
states of the form |4'((7^, 1,1,1)) as in Lemma 3 [case 
(iii)] (see also Table I third row). Let us assume for now 
that none of the state parameters is vanishing, i.e. 7/ 7^ 0 
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and Q 7^ 0, i G {1,2,3}. The necessary and sufficient 
conditions for possible LOCC transformations of these 
states already show, that the corresponding volumes are 
always three-dimensional, as one can vary all three 7} pa¬ 
rameters in accordance with the three inequalities given 
in Lemma 4 [case (iiia)]. Hence, all states |4)(/i^, 1,1,1)) 
are in the accessible set of |'I'(g^,Il,l,l)) if the parame¬ 
ters fulfill 


s.t. |C^| <1/2 and H-r—rfc >0 

for {*,j,fc} = {1,2,3}}, (44) 


where r; = 'yf/Ci- Note that in case with 

w G {x, y, z} additional LOCC transformations are possi¬ 
ble, implying that more states are accessible (see below). 
For that reason we consider first the case for 

any w G {x,y,z}. With the help of cond. (44) one can 
compute the accessible volume of any given state of the 
form |4'(5^, 1,1,1)). It is, however, not easy to give a 
closed expression of the accessible volume for arbitrary 
parameters of g^, as it is difficult to get the limits of 
the integral that gives the accessible volume from cond. 
(44). However, given g^ this volume can be easily com¬ 
puted. Note that for states with 7/ = 0 for exactly one 
k G {1, 2, 3} also states with (1=0 can be reached, see 
Lemma 4 [case (iiib)]. The part of the accessible volume 
containing only states with /■/ = 0 is two dimensional (as 
(I = 0) and hence, of measure zero, as long as the three- 
dimensional accessible volume exists. However, there are 
states |'I'(g^, 1,1,1)) with 7/ = 0 for which cond. (44) 
can not be fulhlled, meaning that they can only reach 
states with (I = 0. Thus, for these states the accessible 
volume is only two-dimensional as explained before (see 
also below for the source volume of states with (I =0). 

Let us now determine the source volume for states of 
the form |4'(g^, 1,1,1)) in terms of the state parameters. 
We consider first the case where none of the coefficients of 
g^ in the Pauli basis vanish. For all states |4)(/i^, 1,1,1)) 
in the source set of |'I'(5^, 1,1,1)) again the three in¬ 
equalities I + Ti — Tj — Tk > 0 for = {1,2,3} 

have to be fulfilled, where here ri = Q/^} (Lemma 4 
[case (iiia)]). These inequalities also imply the necessary 
conditions [y/j > |C/| for I G {1,2,3}. Thus, it can be 
easily seen that the normalization condition on the pa¬ 
rameter vector, i.e. |C^| < 1/2, is automatically fulfilled, 
as [7^ I < 1/2 holds. We can compute the source volume 
of |4'(g^, 1,1,1)) by integrating over the valid parame¬ 
ter ranges for the given by the three inequalities and 


obtain 


H,(|vl/(gi,l,l,l))) = 


^CCi.C^Ca') 

^2 rl+ri—r2 


Pi pr2 pi+ri-r2 

/ / nh2 \l3\ dr3dridr2 

Jo^O —l+ri+r2 

p 1 pl-ri+r2 

/ / nh2\l3\dr3dridr2 

JoJroJ — l+r\ +r2 


f0Jr2'J — 1+ri +r2 

= H\- 


(45) 


Note that out of mathematical convenience we made a 
change of variables in the integral above. Hence, the 
source entanglement is with = l/(36-\/3) given 

by .B«(|4'(5M,1,1))) = 1 - 24V37I7II73II. In Fig. 
7 we plot the source and accessible volume of a state 
|4'(5^, 1,1,1)) of the above form. Both volumes lie in¬ 
side the quadrant of a sphere, which is a part of the 
accessible volume of the seed state, as explained below. 



FIG. 7: (color online) Source (red polyhedron) and accessi¬ 
ble (blue) sets, Ms and Ma (see Eqs. (44) and (45) for the 
volumes), of a 4-qubit accessible state |'F((;^, 1,1,1)) with 
parameters 7} = 0.23, 7} = 0.13 and 7} = 0.15. The vol¬ 
ume of the quadrant of a sphere corresponds to a part of the 
accessible volume of the seed state \'^)^eed- 

Furthermore, the source volume of states of the 
form |4'(5^,1,1,1)]) with vanishing parameters can be 
easily computed with the help of the conditions in 
Lemma 4 [case (iiib) and case (hie)]. If there is one 
parameter vanishing, i.e. without loss of generality 
7} = 0, the two-dimensional source volume is given by 
14(|4'(g^, 1,1,1))) = with the source entanglement 
equal to £’s(|'I'(g^, 1,1,1))) = 1 — Moreover, 

the source volume of a state with two vanishing param¬ 
eters, i.e. without loss of generality 7i = 72 = 0, is 
one-dimensional and reads 14(|'1>((;}, 1,1,1))) = y^ with 
Esil'Sigll,!,!))) = 1-271. 

Note also that there exist states of the form 
|4'((7i, 1,1,1)) for which the accessible volume consists 
of two sets that are not connected (in the used parame¬ 
terization). An example of such a state is given in Fig. 

8 . 
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FIG. 8 : (color online) Source (red polyhedron) and acces¬ 
sible (blue) sets, Mg and Ma, of a 4-qubit accessible state 
|'&((;^, 1 , 1 , 1 )) with parameters — 0.09, 72 = 0.1 and 
73 = 0.08. Interestingly, the accessible set contains two sets 
that are not connected. 


Now we consider the accessible volume of states with 
7 ^ 0 (■"^6 choose again w = x without loss of gen¬ 
erality) and the seed state (see Table I fourth and last 
row). As for the seed states additional LOGO transfor¬ 
mations are possible, we have to treat them separately. 
Let us start with the accessible volume of states of the 
form |l['( 5 ,),, 1, 1 ,1)), with g], 9 ^ 1. We have to split the 
computation of the accessible volume of 
into two parts, as these states can be transformed into 
states of different forms (Table I fourth row). First of 
all according to Lemma 4 [case (iiia)] these states can 
be transformed via LOCC into all states of the form 
|<l)(/i^, 1 , 1 , 1 )) fulfilling the two inequalities 1 ± ri > 0 
[ 35 ] . As we can choose 7i, Ci > 0 the conditions above 
are equivalent to 7i < Ci ■ Furthermore, the two param¬ 
eters (2 ) C3 can be changed arbitrarily as long as they 
obey the condition on the norm, i.e. |C^| < 1/2- Hence, 
the first part of the accessible volume of |dl(g^, 1,1,1)) 
is given by 


r 


1/2.V1/4-(7()" /■V1/4-(7))"-(C3L = 


• 1/2 j-TT gyJTJX^Z 


dC,\dC,\dCX 


= 7r/2 


/O 40 
1 

12 


l\ 


rdrdOdz 

{llf 


(46) 


where we converted the integral into cylindrical co¬ 
ordinates. Furthermore, the state |lE'((/,),, 1,1,1)) can 
be transformed via LOCC into all states of the 
form |$(h;^, hy1,1)) (up to permutations), with 
= 1/21 + Cl see Lemma 4 [case 

(ib)]. The transformation is done by first converting 
|'I'((/^, 1,1,1)) into \^)seed ll^cn converting 

this state into the final state |4>(h;[., 1,1)). Thus, the 

necessary and sufficient conditions on the parameters are 


7 i < Cl (Lemma 4 [case (ib)]) and C 21 C 3 fulfilling the 
normalization condition < 1/2 (with C/ = 0). This 
leads to the second part of the accessible volume given 
by 

/■I /2 /■v'l/4-(C|)^ /■ 1/2 

K.2(|'h(5^,i,i,i)))= / / / dCi'dClrfCl 

J-1/2 Jo 

= ^^(1 - 27 /). (47) 

This part of the accessible volume of |iE'( 5 ,],, 1 ,1,1)) is 
depicted in Fig. 9 for a random state. Furthermore, 
all states given by permutations of the local operators 
of (e.g. |$(/i;^,l,/i 2 _^,l))) are also in 

the accessible set of I'll( 5 ^, 1 , 1 , 1 )) and thus, we have to 
count I 4_2 three times. Therefore, the total accessible 
volume is equal to 

K(|'Ii(gi,l,l,l))) = K.i+3K.2 

= 1/4871(11+ 87iM(7|)'-3]). (48) 

The accessible entanglement is then given by 
Ea{\^{gl 1,1,1))) = 1 + 8 /ll 7 i[( 7 i )2 - 3] as 

ysup _ Furthermore, we have to compute the 

source volume of a state of the form |'l'((/^, 1 , 1 )) 

(Lemma 4 [case (ib)]), as this state can be reached 
on the one hand by states of the form |<f’(/i,]., 1 , 1 , 1 )) 
with C/ < 7 / as we have just seen (Table I fourth row). 
Hence, the first part of the one-dimensional source 
volume of \'^{gl,gl^,„l,l)) is given by 

Vg,^{\'\'[gl,gl,,l,l)))= r\cl=ll (49) 

Jo 

On the other hand according to Lemma 4 [case (ib)] all 
states of the form |'&( 5 ;/, hy .,, 1, 1 )) are also in the source 
set of \'^igl,gl^z, 1,1)) (see Table I first row). Thus, the 
second part of the source volume of \^{g],,gy^^,l,l)) is 
given in Eq. (37), where one simply has to replace 72 ( 3 ) 
with 7 |( 3 ) • The total source volume is obtained by taking 
the sum of both parts and reads 

ys{\'^{gi,gl,^,t,i))) = 7 / + Y^(7|)^ + (7|)^- (50) 

As 1 //“P = 1 the source entanglement is given by 

Esi\'^igl,gl,zAA))) = 1 - ysi\'^igl,gl,zAA)))- It is 
easy to see that the accessible volume of states of the 
form \ ^[g],,gy^^,l,l)) is given by Eq. (35) (see Lemma 
4 [case (ia)]), where one again has to replace 7 / with 7 / 
for i € {1, 2, 3} and set 7 ^ = 0. 

Note that the seed state has the maximum ac¬ 

cessible volume, as we will show in the following. The 
seed state can be transformed up to permutations into 
either states of the form |<l>(/i^, 1 , 1 , 1 )) with arbitrary 
h^, see Lemma 4 [case (hi)], or all states of the form 
|4’(/i(„,Ii|.„,l,l)) with {u,v, w} = {x,y,z}, see Lemma 
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4 [case (ib)] (Table I last row). Hence, the three di¬ 
mensional accessible volume of the seed state is equal 
to four times 1/4 of the volume of a sphere with radius 
1/2 (|^^| < 1/2 and > 0, Ci > 0) plus 36 times 
{{u,v,w} = {x,y,z} and including all party permuta¬ 
tions) 1/2 of the volume of a cylinder with radius 1/2 
and height 1/2, see Fig. 9 (\/(Cu)^ + (C«)^ < 1/2 and 
0 < Ci < 1/2), i.e. 

y^mseed) = V6 + 9/4^ = 29/127r. (51) 

Thus, any seed state obtains the maximum value for the 
accessible entanglement, i.e. Ea{\'^)= 1. 



FIG. 9: (color online) A part of the accessible set (blue), 
Ma (see Eq. (47) for the volume), of a four-qubit state 
|4'((;/, 1,1,1)) with 7 i = 0.2. The volume is equal to 1/2 the 
volume of a cylinder with height ( 1/2 — 71 ) and radius 1 / 2 . 
The red line is the source set of |4'((;^, 1,1,1)) (the source 
volume is one-dimensional). The volume of the larger set cor¬ 
responds to a part of the cylinder belonging to the accessible 
volume of the seed state, for which 7 } = 0 . 


Summarizing the results on the accessible and source 
entanglement of generic four-qubit states we highlight 
again the three-dimensional accessible volume of the seed 
states, which is the maximum accessible volume of all 
generic states. This can already be seen in Table I, where 
in the last row for a conversion from the seed state to 
other states no additional necessary and sufficient con¬ 
ditions apart from the condition on the norm of the pa¬ 
rameter vectors have to be fulfilled. Hence, the seed 
states can reach the most other states deterministically 
via LOCC. Furthermore, when considering transforma¬ 
tions between states as in the first row of Table I the 
states in the MES, i.e. have a two- 

dimensional accessible volume and they maximize obvi¬ 
ously the source entanglement. Whereas, states that are 
not in the MES have in this case only a one-dimensional 
accessible volume and are therefore, infinitely less pow¬ 
erful than states in the MES. 


V. CONCLUSIONS 

In [14] we have introduced two classes of entangle¬ 
ment measures for general multipartite quantum states: 
the accessible entanglement and the source entanglement. 
These measures have a clear operational meaning related 
to the relative usefulness of the state under single-copy 
deterministic LOCC manipulation. Moreover, whenever 
these transformations can be characterized, the measures 
can be computed. In [14] we showed this in the case of 
three-qubit pure states and provided explicit formulae in 
these cases for both Ea and Eg. In this article we have 
further demonstrated the applicability of our approach 
by analyzing the case of bipartite pure states of arbitrary 
dimension and of generic four-qubit pure states. In the 
first case, using tools of convex polytopes, we provided a 
closed expression for the source entanglement of an arbi¬ 
trary pure bipartite state of Schmidt rank d considering 
all possible dimensions on which Mg can be supported, 
{Eg}k>d- The accessible volume turns out to be more 
complicated; however, there exist algorithms that allow 
to compute all measures {E^}k<d- Moreover, in the case 
of bipartite states of low dimension such as d = 2,3 we 
have also obtained analytic expressions for the accessible 
entanglement. In the second case of four-qubit generic 
states, we have completed the analysis of [9] characteriz¬ 
ing all possible LOCC conversions inside this class. Using 
this we have provided explicit formulae for Ea and Eg for 
all these states. 

The results presented here show the versatility and 
applicability of the measures we introduced in [14], 
having covered the most relevant classes of pure few- 
body entangled states. We hope that having measures 
that are both computable and with a clear operational 
meaning will help to understand better the properties 
and potential applications of multipartite quantum 
systems. For the future, our research opens several 
directions for further investigation. First, it seems 
desirable to analyze in detail more many-body cases 
with the aim of obtaining closed expressions for our 
measures (or efficient algorithms for their computation). 
In this context it would also be interesting to exploit the 
results of [27] to try to evaluate our measures for mixed 
states. Moreover, we would like to connect our measures 
with applications. An interesting option is to study 
the role of our measures as figures of merit for known 
quantum information protocols. This could then lead 
to the identification of the most relevant multipartite 
states and maybe allow to devise new applications of 
multipartite entanglement. Also, we hope that our 
measures can bring new insights to understand the 
interplay between entanglement and many-body physics 
in, for example, phase transitions. Apart from that, the 
investigation of LOCC transformations among several 
copies of a given state and e-LOCC transformations [36] 
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in this context will be relevant. 

After completing this work we became aware of Ref. 
[28], where the volume of polytopes of the form of 
in Eq. (13) was calculated. 

This research was funded by the Austrian Science 
Fund (FWF) Grant No. Y535-N16 and the Spanish 
MINECO through grants MTM 2010-21186-C02-02 and 
MTM 2011-26912 and the CAM regional research con¬ 
sortium QUITEMAD-hCM S2013/ICE-2801. 


Appendix A: The source volume of non-generic 
states 

In this appendix we show that the formula for the 
source volume given in Eq. (19), which was derived for 
states with non-degenerate Schmidt coefficients, is valid 
in general. This derivation was based on the fact that 
is a simple polytope for a generic state, such 
that the formula in Eq. (15) was applicable to compute 
its volume. However, {iji) fails to be simple in gen¬ 
eral for degenerate states, i.e. states with at least two 
identical Schmidt coefficients. Consider, e.g., the state 
with Schmidt vector A('!/') = (Ai, A 2 , A 2 , A 4 ) for which 
A 2 = A 3 . This vertex fulfills the four inequalities 

El{\') < Eu{Xm, fee { 1 , 2 , 3} 

out of the set 17, which defines (see Eq. (11)), 

with equality. Hence, X{tjj) is an element of more than 
d — 1 = 3, more precisely of four, facets which in turn 
implies that the three-dimensional polytope is not simple. 
It is also straightforward to show that X{'tjj) has more than 
d—1 = 3 neighbors. Indeed, each of the following vertices 

(A 2 , A 2 , Ai, A 4 ), (Ai, A 4 , A 2 , A 2 ), 

(A 2 , Ai, A 2 , A 4 ), (Ai, A 2 , A 4 , A 2 ), 


is an element of d — 2 = 2 facets that also contain X{ip) 
and is thus a neighbor of X{ip)- We therefore conclude 
that Eq. (15) is not directly applicable to non-generic 
states. 


However, in deriving Eq. (19) one obtains that 


VsW = 

1 y/d -pr |Afc —Afc+ij 

^k-Xk+i 


fc=i 



d-l 


d-1 


n 


cr(fc -I- 1) 
(AI) 


Due to the denominator in the first product of Eq. (AI), 
degenerate states must be treated with care in deriv¬ 
ing a simple expression for the source volume. We show 
that, for any direction (along non-degenerate states) from 
which we approach A('0) within the subset of sorted 
Schmidt vectors, we obtain in the limit Vs{ip) (given by 
Eq. (19)) as the volume of the source set. Stated dif¬ 
ferently, for any non-degenerate, sorted Schmidt vector 
A and any e > 0 we define the non-degenerate, sorted 
Schmidt vector A(A, e) = (1 — e)X{tj;) -|- eX for which it 
holds that 


limH,(A(A,e)) = H,(^), (A2) 


with I4(V’) given by Eq. (19). Here, we used that the 
set of ordered Schmidt vectors is convex and that A(A, e) 
is a non-degenerate Schmidt vector for which Eq. (AI) 
is defined. Let us now show that Eq. (A2) is valid. We 
denote by I the set of all indices i G ( 1 ,..., d} for which 
Xi = Ai+i. Then we obtain the following. 


lim t4(A(A, e)) = 
£->■0 


1 . 

—r ^-r7 hm 

dl {d — 1)! e-s-O 


n 

fc^i 


(1 ~ e)(Afc — Afc+i) -l- e(Afc — Afc+i) 


n 


-{Xk — Afc+i) 


E 


(1 — e)(Afc — Afc+i)-I-e(Afc — A/c+i) e(Afe — A^+i) a^Sd 


1 Vd 


E 


E cr(fc)Afc - ^ 

k=l 


i-1 


d\(d-l)\ ^ d-l 

<y&Sd J)[ cr{k) — a{k + 1) 

k=l 


= Vs{iP). 


/ d _ .d-l' 

(E^(fc)A(A,e)fc-^j 


n cr{k) — (j{k + 1) 
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As e, Afe — Afe+i > 0 the second product gives 1 and the 
limit can be computed easily. We thus proved that for¬ 
mula Eq. (19) for the source volume also holds for de¬ 
generate states and is therefore applicable in general. 

An example of a degenerate state is the maximally en¬ 
tangled state with Schmidt vector IjVd- . Since 
the maximally entangled state cannot be obtained from 
any other state via LOCC its source volume has to be 
zero. Indeed, it holds that 

E^ ^ E- ^ = 0 . 

fc=l “ k=l 

If one inserts this expression into Eq. (19) it is easy to 
see that 14 (</>''’) = 0 . 


Appendix B: Two different approaches to compnte 
the source and accessible volume 

In the main part we used two different ways to calcu¬ 
late the (d—I)-dimensional volume of convex polytopes in 
the space A C R‘^. In the first one, used in the derivation 
of the source volume, we directly determined the volume 
by considering the d-dimensional Schmidt vectors in this 
(d— l)-dimensional subspace oiR^. However, in order to 
determine the accessible volume we considered the pro¬ 
jection of the accessible volume onto the subspace of 
that is spanned by the first (d—I) standard vectors. That 
is, we treated the accessible polytope as a subset of the 
set {(Ai,..., Ad_i)|Ai > A < 1} However, 

in order to use the same measure for both, the source 
and the accessible volume, the volume element in the 
subspace that contains the projection of the latter has to 
be multiplied by the Jacobian of the coordinate transfor¬ 
mation that relates the two methods. As we explain in 
the following, this amounts to multiplying the volume of 
the projection of the accessible set by a constant factor 

of Vd. 

We have to determine how a volume element ex¬ 
pressed in an orthonormal basis (ONB) of ?7 = {A £ 

R‘^ s.t. = 0} related to a volume element in the 

projection oiiR‘^~^, i.e. the subspace on which we project 
the accessible set. It is straightforward to show that the 
following vectors are an ONB of R‘^, 

fio = 

Aifc = ^JL_ (I,..., I, -fc,0,..., 0) for fc G {1,... ,d- 1} 

where in a total of k entries of 1 are followed by 
one entry of —k and zeros. The vectors are 

also an ONB of U. If a vector X G R‘^ has coordinates 
(mq. Ml,..., Ud-i)^ in this ONB, its coordinates in the 
standard basis A = (Ai,..., Ad-i, 1 —given 


by 


d-l 

X = '^Ukflk, (Bl) 

k^O 

where Ug = Ij^fd for a normalized Schmidt vector. It 
is therefore easy to find the d — I Schmidt coefficients 
{Ai} as a function of the coefficients {Mi}. A volume ele¬ 
ment in U is given by Of=i while a volume element 
in the projected space is given by 11^=1 B is well- 
known that these volumes are related by the Jacobian, 

TulZ-ZZi) ’ according to 


d-l 

= 

i=l 


det 


/9(Ai,..., Xd-i) 

\d{ui,.. .,Ud-i) 



(B2) 


Using Eq. (Bl) it is easy to show that 


det 


/ 5(Ai,..., Arf_i) 

V9(Mi,...,Md-l) 


d-l 

n 


\/k^/k + 1 


det 


(I 

-1 

0 

0 



VO 


1 1 ... l\ 

1 1 ... 1 

-2 1 ... 1 

0 -3 ... 1 

0 0 ... ly 


Hence, we have that 0 ^= 1 ^ dui and the 

volume of the projection is v^-times smaller than the 
original volume. However, since this is a constant factor 
and since we rescale the volumes in order to obtain the 
entanglement measures, the two different ways to calcu¬ 
late the volumes of the convex polytopes lead to the same 
result for the entanglement measures. 


Appendix C: Proof of LOCC convertibility 
conditions 

In [9] many necessary and sufficient conditions for 
LOCC transformations of generic four-qubit states were 
obtained (in particular case (i) and case (iii) have been 
considered there). Whereas in [9] only the conditions 
for a state to be reachable (or convertible) have been 
derived. Lemma 4 gives the necessary and sufficient con¬ 
ditions for the existence of a LOCC-protocol. In order to 
prove Lemma 4 we use that LOCC is strictly contained 
in the set of separable (SEP) operations, for which the 
criterion for their existence for state transformations has 
been derived in [26]. Here, we first review these results 
and then use them to prove Lemma 4. 
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Let us denote by S{'i'seed) = {S £ G : = 

l^)seed} of symmetries of the seed state, where 

G denotes the set of local invertible operators. Then for 
g,h £ G the state IT) oc can be transformed 

into the state |T) oc via SEP iff there exists a 

m gM, a set of probabilities and Sk G 5'(Tseed) 

such that [26] 

Y,PkSlHSk=rG. (Cl) 

k 

Here, H = h^h, G = g^g are local operators and r = 
II l'^) IP/II 1^) IP- This criterion is used in [9] to derive the 
following necessary conditions for LOCC convertibility of 
generic four-qubit states. For these states given in Eq. 
(32) there are finitely many symmetries of the seed states 
given by {<7i}^^Q. Hence, Eq. (Cl) implies [9] 

Si{H) = (C2) 

where we used the same notation as in [9], i.e. H = 
Here, £4 is the completely positive map given 
by Eq. (Cl), i.e. £i{H) = J2kPk<^k*Haf^ and £i{p) = 
'^kPk^k^ P^k^ ■ Note that by taking the trace of Eq. (Cl) 
one obtains r = 1. Furthermore, by tracing over the last 
two subsystems in Eq. (C3), one obtains £ 2 {H^ 0 H'^) = 
£i{H^) ® £i{H'^) (and similarly for other subsystems). 
Whereas this equation has been used in [9] to identify 
the reachable states, here, we mainly use the necessary 
condition obtained by tracing over all but one system, 
i.e. 

£i{H^)=G\ (C3) 

which is equal to rfQC^ = 7 ^. We use again the same 
notation as in [9] with Q denoting the Hadamard prod¬ 
uct (i.e. entry-wise multiplication), fj = (771,772,??3)^, 
Vo = J2k=oPk = 1 and 77 * = pg + K - (Pj + Pfe), with 
{i,j, fc} = {1,2,3}. Hence, we use the necessary and 
sufficient condition for SEP transformations in Eq. (Cl) 
to obtain necessary conditions for the LOCC convertibil¬ 
ity, Eq. (C3), and show that they are also sufficient by 
constructing a corresponding LOCC protocol. 

1. Proof of Observation 2 

Due to the fact that the non-isolated generic four- 
qubit states are of the form [9], with g £ 

U^ue{x,y,z}^^’ Gw = {g^ ® gl® gi® gt„gl, ® 
g'^ ® gi® gt,gi ® gl ® g^ ® gt,gi ® gi ® gi ® g'^}, 

we only need to consider transformations among these 
states. First we show that g with g £ Gw can 

only be transformed to a state h £ Gw for the same 
value of w. This can be easily seen by contradiction. 
Suppose that h £ Gv with v ^ w. There always ex¬ 
ists an i € {1,2,3,4} for which Eq. (C3) implies that 


£i{Hl) = Yl,kPk<^kHlcrk = G\^. However this cannot be 
fulfilled ioY V ^ w as <7kHl<7k G span{l, (T„} for any fc, 
whereas G\^ £ span{I, ctu,}. Hence, the transformation 
above is only possible if v = w. For the remainder of 
the proofs we will therefore set w = v = a; in order to 
simplify notations. The other cases, i.e. v = y,z can 
be straightforwardly obtained from that. A similar argu¬ 
ment as above shows that the position of the operators 
which do not only have a component in x-direction must 
coincide. That is, e.g. a transformation of the form 

gl®g^®gl®gt \'^)seed h^® hi® hi® hi 

(C4) 

is not possible. This completes the proof of Observation 
2 . □ 
Hence, due to Observation 2, the only possible transfor¬ 
mations that have to be investigated for LOCC are (up 
to permutations of the parties) 

g^®gl®gl®gt\^)seed h^®hl®hl®hl\'^)^^^^. 

(C5) 

2. Proof of Lemma 4 

Note again that Lemma 4 [case (ia)] has been proven 
already in [9]. Out of convenience we proof Lemma 4 
[case (ib)] at the end of this appendix, as we use some 
of the results obtained in the proofs of the other cases 
there. Hence, we go directly to Lemma 4 [case (ii)], where 
transformations of the form (up to permutations) 

gl®gl®l®l hl®hl®l®l 

(C 6 ) 

are considered. Obviously, the proof works analogously 
if we choose gl instead of g^. It follows from Eq. (C3), 
i.e. the condition £i{Hl) = Gj, that 7 ^ 2(2 = I 2 
thus, 7 } can only be increased, i.e. 7 } < Cl- We 
determine a simple two-outcome POVM that realizes 
the transformation. The POVM applied by the first 
party is given by {v^h),(gl)”\ - phlcr^{gl)-^} with 

p £ [0,1]. This is indeed a POVM if (2p — l)Cl = 7 } 
and can be implemented by LOCC as for the second out¬ 
come all other parties simply apply the LU <7x- After 
applying this POVM the state |T(hy,g}, 1,1)) is ob¬ 
tained. The second transformation works analogously. 
Again, the parameter 7 ^ can only be increased by LOCC, 
i.e. < Cl- The second party applies the POVM 
{VP^l{9l)~^, V'^-phlcTyigl)-^}, which can also be im¬ 
plemented by LOCC and thus, Lemma 4 [case (ii)] is 
proven. 

For transformations of the form (Lemma 4 [case (iiia)]) 

0 1 0 1 0 1 010101 , (C7) 

the necessary and sufficient condition for the existence of 
such a separable transformations can be easily seen to be 
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given by (see Eq. (C3)) 

G^=£l{H^) = l/21 + Y,Clv^<^^, (C 8 ) 

i 

where as defined before t]o = Ylt=oPk = 1 iji = 
Po +Pi - {Pj +Pk), with = {1,2,3} and pi > 0 

for any i. Due to the condition above we have for C,} ^ 0 
for all i, Pi = 1/4(1 + — rj — rk) for = jl/Ch where 

{i,j, k} = {1, 2,3} and po = 1/4(1 + ri + r 2 + r^). Hence 
the necessary and sufficient conditions for the existence 
of a separable transformation are that all these prob¬ 
abilities are non-negative, as stated in Lemma 4 [case 
(iiia)]. Note that any of these separable transformations 
can be realized via LOCC as the POVM {Mi}f^^ [9] 
with Mi = y/pih^(Ji{g^)~^ 0 (tTi)®^ can be implemented 
locally. Hence, if Q ^ 0 for all i the necessary and suffi¬ 
cient conditions for the existence of the LOCC transfor¬ 
mation are 1 + ri — Vj — rk >0 for any choice of 
such that {i, j, fc} = (1, 2, 3}. 

Note that we have different necessary and sufficient 
conditions for transformations as in Eq. (C7) if the 
state |<I)(/i^, 1 , 1 , 1 )) has vanishing (// parameters as 
then ri above is not defined. Eor states with one 
vanishing parameter (Lemma 4 [case (iiib)], i.e. without 
loss of generality Ci = 0, Eq. (C 8 ) implies that also 
the corresponding 7 / parameter has to be equal to 
zero, i.e. 7 } = 0. The other two parameters have to 
fulfill ?72(3)C2(3) = 7 ^( 3 ) with 77 , = Pq+P i - (pj + Pk), 
{i,j,k} = {1,2,3}. As r ]2 and 73 are linearly in¬ 
dependent and the two parameters 7 }, 7 ^ can be 
chosen nonnegative in the standard form (Eq. (32)) 
the state |<i)(/i^, 1 , 1 , 1 )) can be reached from all states 
|'I'(( 7 ^, 1 , 1 , 1 )) with 7 } < and 7 } < Cs- A similar 
condition holds for states |$(/i^, 1 , 1 , 1 )) with two 
vanishing parameters, i.e. without loss of generality 
Cl = C2 = 0- Then, from Eq. (C 8 ) it follows again 
that the corresponding 7 / parameters are equal to 
zero, i.e. 7i = 7 } = 0. Furthermore, the condition 
VsCs = 73 tias to be fulfilled, thus, the only necessary 
and sufficient conditions for transforming |'I'(( 7 }, 1 , 1 , 1 )) 
into |$(/i}, 1 , 1 , 1 )) is 7 { < C{. 

Note that the necessary conditions for reaching the 
state |$(/i},, 1,1)) via LOCC (Lemma 4 [case 

(ib)]) can be obtained by considering the equation 
£2{Hx®Hy .,) = £i{Hl)®£.i{Hy^^) (see [9]). From there 
it follows that 771172 — 773 = 0 and 771773 — 772 = 0. Hence, 
these conditions are on the one hand fulfilled if (771)^ = 1, 
which leads to the same necessary and sufficient con¬ 
ditions as in Lemma 4 [case (ia)]. The states fulfilling 
these conditions are of the form |'I'(/i}., 1 , 1 )) and 
a POVM allowing for this transformation via LOCC is 
given in [9]. On the other hand the necessary conditions 
are also fulfilled if 772 = 773 = 0. By taking into account 
Eq. (C 3 ) this implies that 771^} = 7} and 7! = 73 = 0. 
Thus, |<i)(/i}.,1, 1 )) can also be reached by states 


of the form 147 ( 5 },, 1,1,1)) with 7} < C}. The POVM 
realizing this transformation via LOCC are already 
given in [9] and in the proof of Lemma 4 [case (iii)] 
above, as we simply first transform |47(5}, 1,1,1)) into 
|4>(/i}, 1,1,1)) and then the latter state is converted 
into the final state |$(/i}, l-i I))- Cl 
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[33] This can be easily seen as follows. We take an ar¬ 

bitrary vertex \{tp)a of the source set. This vertex 
fulfills the inequalities {X{’ip)a) < Ek{X{ip)),k G 

{1,..., d — 1}, as an equality. Based on our previous re¬ 
sults we therefore conclude that the vertices of the form 
P<yPTi,i+iPa-T^P<yX{’il)) = PaPr^ ^^^X{'<|J) = X{'tl))a,i are the 
only neighboring vertices of Xi^l))^ for k G {1 ,..., d — 1}. 

[34] Indeed, it holds for any matrix A G that ]det(A)[ = 
]det(l©A)], where the left-hand-side is the volume of 
a fc-dimensional parallelotope defined by the column- 
vectors of A, and the right-hand-side is the volume of a 
{k -b l)-dimensional polytope that has the previous poly¬ 
tope as base to which the first vertex is orthogonal. 

[35] Note that in this case r 2 = ra = 0. 

[36] With e-LOGC transformations we refer to LOCG trans¬ 
formations in which a small deviation in the target state 
is allowed with respect to some suitably chosen distin- 
guishability measure. 


